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STEADY-STATE ANALYSIS OF A BRAYTON SPACE-POWER SYSTEM 

by john L. Klann 
Lewis Research Center 


SUMMARY 

An experimental Brayton space-power system was approximated in a digital com- 
puter study. Parametric and simulated-fixed- system operations to be conducted in the 
experimental system ground tests were examined. The study shows trends in system 
conversion efficiency for both design and off-design operating conditions. The analysis 
is independent of operational limits. However, the present limitations of the existing 
hardware are discussed. 

With the design working-gas mixture of helium and xenon, and with temperatures at 
their design values, estimated levels of system conversion efficiency increased from 
about 0. 14 at 2 kilowatts of output to about 0. 24 in the 10- to 15-kilowatt range. Changes 
in conversion efficiency with off-design turbine- and compressor-inlet temperatures 
were calculated. Each inlet temperature was individually varied while the other was 
maintained at its design value. System conversion efficiency was insensitive to ±20 per- 
cent changes in the relative amount of helium in the gas mixture. However, the use of 
krypton in place of the design gas mixture resulted in up to a 0. 03 decrease in conver- 
sion efficiency. 

Study of simulated radioisotope-power- system operations showed that bleeding- 
working gas would be an effective means to maintain high system power outputs during 
the decay of a short-lived isotope. 

Future power system gains could be achieved by modifying the compressor-diffuser 
vane angle for peak compressor performance. With the modified diffuser, the analysis 
showed conversion efficiency gains of 0. 01 at the 2-kilowatt power level and of 0. 02 in 
the 10- to 15-kilowatt range. 


INTRODUCTION 

Lewis Research Center has been advancing closed-Brayton-cycle technology for 
space electric-power-generation systems since about 1963. Current efforts are cen- 



tered. on a Brayton power system intended for the output power range from 2 to 10 kilo- 
watts. The system has a design- life goal of 5 years. Its design working gas is a mix- 
ture of helium and xenon at the molecular weight of krypton, 83. 8. Within the closed 
system, a single- shaft assembly of small rotating machinery, supported by gas bearings, 
operates at a controlled speed of 36 000 rpm (3770 rad/sec) and produces electric power 
at 1200 hertz. 

The cycle operating parameters for this system were selected in reference 1 on the 

basis of minimizing total system weight with a radioisotope heat source. The weight of 

the radioisotope heat source and its shield was dominant. To keep the source size and 

weight low, the Brayton-cycle operating parameters were chosen to favor efficiency at 
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the expense of increased specific radiator area (ft /kW) for cycle waste-heat rejection. 
Although the system design conditions were selected for a radioisotope heat source, the 
system could be mated to either a nuclear reactor or a solar-heat collector. 

Lewis is currently ground testing this 2- to 10-kilowatt power system. Hardware 
components for the system were built under contracts with the aerospace industry. 
Reference 2 describes the components and the intended system ground-test configura- 
tions. In addition to the system tests, several of the major components and subsystems 
are being investigated separately. References 3 to 6 present some of these results. 

The first system ground test is using a radiant electric heat source for parametric study 
of system performance and to simulate both solar- and radioisotope- power- system op- 
eration. The ground testing also simulates system space-heat rejection with the use of 
test-facility support equipment. 

The electrically heated ground-test configuration was approximated in a digital 
computer study of its design and off-design operation. Initial results of this study are 
presented herein. The power system configuration was synthesized. Where measured 
component performance was not available, design predictions were used. Approximate 
heat-transfer and pressure-drop relations were developed, fitted to the component per- 
formance, and used in the analysis. 

Prior to presenting the analytical results of the study, the present power system is 
described, and present limitations of the existing hardware are discussed. However, 
the analysis itself was made independent of any system limits, and output power ranges 
from 2 to about 15 kilowatts were examined. 

The study was made mainly to show the trends in system conversion efficiency to be 
expected from the ground-test configuration. Ranges in turbine- inlet temperature, 
compressor-inlet temperature, and the relative helium mass in the working-gas mix- 
ture were examined. The use of krypton in place of the design gas mixture was studied. 
Calculations were made for simulated solar- and radioisotope- power-system operation. 

A promising future change to the compressor-diffuser vane angle was studied for its ef- 
fect on system operating conditions. 
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POWER SYSTEM DESCRIPTION 


Test Configuration 

Figure 1 is a schematic diagram of the system test configuration. The station 
numbers and all symbols are defined in appendix A. 

Heat source. - The heat source consists of two planar banks of quartz lamps with a 
central, U-tube heat exchanger. This heater can be operated in either a constant- 
turbine- inlet-temperature mode for parametric studies and simulation of a solar-heat- 
collection system, or in a constant-gas-thermal-input mode to simulate a radioisotope 
heat source. 

Rotating machinery . - The rotating machinery (turbine, alternator, and compres- 
sor) forms a single- shaft assembly that is supported by gas bearings. The machinery 
is designed for constant- speed, steady-state, operation at 36 000 rpm (3770 rad/sec). 
The turbine and compressor are single-stage, radial-flow machines, and the alternator 
rotor is a solid, brushless, four-pole machine of the Lundell type. A speed control 
and a parasitic, resistive load are used to maintain a constant output frequency of 
1200 hertz and to absorb unused power. A regulator maintains the three-phase voltage 
at 120 volts, line-to- neutral, or 208 volts, line-to-line. 

Recuperator and gas-to-liquid heat exchangers . - In addition to the source heat 
exchanger, there is a recuperator and a gas-to-liquid, heat-sink, heat exchanger (or 
waste heat exchanger) in the closed working-gas loop. The recuperator is a single- 
pass, counterflow heat exchanger with plate-fin surfaces. The waste heat exchanger 
has a cross- counterflow arrangement with plate-fin surfaces for both flows. The liquid 
coolant makes eight passes across the straight-through gas flow. There are two sep- 
arate and completely independent liquid circuits, only one of which is to be used at any 
given time. 

Electrical subsystem . - The electrical subsystem includes the speed control, the 
voltage regulator, the parasitic load, a direct-current power supply, two direct-current 
to 400-hertz inverters, and a signal conditioner. During steady- state power system 
operation, the direct-current power supply rectifies part of the alternator output to 
provide internal needs at ±28 volts. Direct- current power is used by the engine con- 
trols, the signal conditioner, and by either of the two inverters. The function of the in- 
verter is to provide 400-hertz power to a motor driving a liquid pump in the heat- 
rejection subsystem. 

Heat-rejection subsystem . - The heat-rejection subsystem contains dual coolant 
loops, each with a separate motor-driven pump. The two completely independent loops 
provide total cooling redundancy. The coolant is dimethyl polysiloxane, 
(CHgjgSifoS^CHgjglgOSifCHgjg, a silicone liquid with a room-temperature viscosity of 
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Parasitic- load 



Facility refrigerant 


Figure 1. - Schematic diagram of Brayton-power-system test configuration. 
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2 centistokes (2x10” m /sec). In the test arrangement, the liquid cools three compo- 
nents in parallel (see fig. 1). They are the waste heat exchanger, the alternator, and a 
series of four flat-plate heat exchangers (or "cold plates"). Electrical subsystem com- 
ponents are mounted on these cold plates. The three flows mix at the exchanger outlets 
and pass through one of two test- support heat exchangers. After being cooled, the liquid 
returns to the pump for recirculation. All of the flow from each pump goes through the 
pump-motor cavities to cool the motor and lubricate its bearings. 

The test- support heat exchangers and facility cooling supply, along with a dynamic 
controller, provide variable waste-heat-rejection capacity. The controller can be used 
to provide facility cooling in either of two modes. One mode results in a constant 
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compressor-inlet gas temperature. This mode will be used for the study of system 
parametric operation. The other mode results in waste-heat- exchanger, liquid- coolant 
temperatures that would occur in space with a fixed-area radiator. This mode will be 
used for the study of simulated solar- and radioisotope-power- system operation. 

Gas management subsystem. - The gas management subsystem supplies and controls 
the power system working gas. It is used during power system starts and stops and for 
gas inventory control. 


System Design Conditions 

Table I shows the parametric values used in the power system component designs. 

All of these values, except for compressor specific speed, were taken from reference 1. 

The design compressor specific speed of 0. 698 was used, rather than the value of 0. 775 

given in reference 1. This resulted in increased working-gas supply pressure levels 

available to the machinery bearings for self-acting operation. 

Any power system must, on demand, be capable of producing power from zero to its 

rated output. The rated output of a closed-cycle Brayton system may also be adjusted 

by changing the working-gas inventory. Gas inventory, pressure level, and mass-flow 

rate are proportional to each other. With a constant shaft speed, the shaft work and 

alternator output are proportional to the gas mass-flow rate. Hence, by setting the gas 

pressure level the rated output of the system may be adjusted. 

All components for this power system were designed to operate over the system 

output power range from 2 to at least 10 kilowatts. The design compressor-outlet pres- 

2 

sure range, at a constant pressure ratio, was from about 14 to 45 psia (97 to 310 kN/m 
abs). These pressures were calculated to result in from 2 to 10 kilowatts of system 
output power. The design effectiveness parameters (table I) and pressure drops for the 
recuperator and waste heat exchanger were used at the estimated 10-kilowatt, gas-loop 
conditions. However, in contrast to the rotating- machinery maximum design pressure 


TABLE I. - DESIGN PARAMETER LIST 


Working gas 

Helium and xenon mixture 

Working-gas molecular weight 

83. 8 

Turbine- inlet temperature, °R (K) 

2060 (1144) 

Compressor-inlet temperature, °R (K) 

540 (300) 

Shaft rotational speed, rpm (rad/sec) 

36 000 (3770) 

Compressor pressure ratio, r^. 

1. 90 

Compressor specific speed, N g ^ 

0.698 

Cycle-loss pressure ratio, 

0. 92 

Recuperator effectiveness, E 

0. 95 

Waste-heat-exchanger effectiveness, E w 

0. 95 

Waste-heat-exchanger capacity-rate ratio, C w 

0. 87 
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of 45 psia, the maximum heat- exchange!' design pressure was 56 psia (386 kN/m abs) 
at the recuperator cold- side inlet. 

The rotating machinery, recuperator, and waste heat exchanger were designed for 
continuous, high-temperature operation over a 5-year period. 


Present Hardware Limits 

For system operation at design temperatures and at constant machinery rotational 
speed, the limits on existing components are functions of gas pressure level and time. 
However, these operational limits are not yet well defined. The low-pressure (or low- 
power) operational limit would be due to insufficient bearing- cavity pressure for self- 
acting operation. This limit is not time dependent. For now, it is assumed that the 
low-pressure limit is at the minimum machinery design operating pressure level, which 
corresponds to a compressor-outlet pressure of about 14 psia (97 kN/m abs) and a 
power output of 2 kilowatts. 

The high-pressure (or high-power) operational limit would probably be due to 
(1) limited capability for alternator cooling, (2) limited thrust-bearing capacity, or 
(3) a structural limit in the heat exchangers. Each of these potential limits may occur, 
in turn, at particular values of high operating gas pressure levels. Both alternator- 
cooling limits and heat- exchanger- structural limits are time dependent. Examination 
of design information and preliminary component test results indicates that for the sys- 
tem design- life goal of 5 years, a limit on alternator heating would occur first. 

Based on unpublished statistical data for alternators with similar windings and 
potting compounds, a 5-year life corresponds to a winding temperature of about 490 K. 
Alternator tests in room-temperature air (ref. 4) showed that the winding hot-spot tem- 
perature increased as the load power factor was reduced from 1. 0 to 0. 75. At a gross 
alternator output of 13. 4 kilowatts with a power factor of 0. 75, the hot-spot temperature 
did not reach more than about 460 K. However, the machinery design environment was 
not duplicated. 

Preliminary, unpublished results from a test of a single set of the rotating machin- 
ery at design conditions showed that the alternator hot-spot temperature reached 493 K 

at 11.4 kilowatts of gross alternator power with a load power factor of about 1. 0. The 
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corresponding compressor-outlet pressure, about 38 psia (262 kN/m abs), was below 
both the machinery and heat- exchanger maximum design pressures. The machinery test 
was continued at design cooling flow conditions until, at a gross alternator power output 
of about 12. 7 kilowatts, the alternator hot-spot temperature reached 513 K. Estimated 
winding potting compound life at 513 K would be about 1 year. With additional cooling, 



a gross alternator output of 15 kilowatts with a hot-spot temperature of about 530 K was 
obtained in the machinery test. 

A theoretical check on the load capacity and film thicknesses of the thrust bearing 

showed its operation to be practical but marginal at a compressor-outlet pressure of 

56 psia (386 kN/m abs), the maximum heat- exchanger design pressure. Preliminary 

test results, however, have indicated lower-than-design thrust loads, so thrust-bearing 

operation above 56 psia (386 kN/m abs) may be practical. Analytically, the power 
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system was not explored above 65 psia (449 kN/m abs). 

METHOD OF ANALYSIS 

A digital computer program was developed to study the steady- state operation of 
the electrically heated ground-test power system. The configuration was synthesized. 
Data from individual rotating-machinery component tests were used. Performance 
maps for the turbine were from Nusbaum and Kofskey's results (ref. 3). Alternator ef- 
ficiency was obtained from Repas and Edkin's work (ref. 4). The compressor perform- 
ance maps were obtained from unpublished results of tests conducted at Lewis. Mathe- 
matical models were developed for the heat-transfer components. These models were 
fitted to either design predictions or to design predictions modified by acceptance test 
results. Details of the computer methods are presented in appendix B. Component 
mathematical models and their use are presented in appendix C. 

The analysis was made independent of any system limits over output power ranges 
from 2 to about 15 kilowatts. The entire analysis was made only at the design machin- 
ery physical speed of 36 000 rpm (3770 rad/sec). Perfect insulation was assumed, so 
there was no heat leakage from the ducting, heat exchangers, or rotating-machinery 
housing. Component performance, used as the basis for the analysis, did not include 
the effects of flow maldistributions due to system ducting. Hence, uniform gas-flow 
distributions were implied in the analysis. 


Gas Loop and Output Power Calculations 

For each computer calculation, the closed gas loop was constrained to satisfy 
values of turbine-inlet temperature, compressor- inlet temperature, and working-gas 
inventory. With these constraints a series of iterations was needed (appendix B). Each 
calculation located those points on the turbine and compressor performance maps which 
satisfy the closed-loop pressure drops. During the course of finding the operating con- 
ditions, both the required ideal thermal input to the gas in the source heat exchanger 
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and the gross shaft power were determined. The gross shaft power includes a penalty 
for shaft-seal losses (appendix C) but not bearing friction (eq. (C34)) and alternator 
windage losses (eq. (C32)), which are subtracted from the gross shaft power. The 
gross alternator output was calculated (appendix C) for an electromagnetic efficiency 
corresponding to a load power factor of 0. 75. Internal system electric power needs 
were approximated (appendix C and eq. (C100)) and subtracted from the gross alternator 
output to obtain the gross system output power P gy G> Gross system output power is 
defined as the sum of the three-phase, 1200-hertz power available at the parasitic- and 
user-load bus (see fig. 1). 


Liquid-Coolant-Loop Calcu lations 

The heat loads from the alternator cooling and the electronic- component cold plates 
(fig. 1) were not included in the analysis. These heat loads require lower coolant- outlet 
temperatures than that from the waste heat exchanger (ref. 2). In a space version of the 
power system, a user might reject these heat loads and further power- conditioning heat 
loads in a low-temperature radiator. The alternator and cold-plate cooling needs do not 
affect the power output calculations and are beyond the scope of this analysis. 

For the study of parametric system operation, the compressor- inlet temperature, 
turbine- inlet temperature, and gas inventory were assigned for each calculation as in- 
dependent parameters. After gas loop conditions were satisfied, the coolant mass-flow 
rate was calculated over a range of values (appendix B). At each value of coolant mass- 
flow rate, the particular values of waste- heat- exchanger coolant-inlet and outlet tem- 
peratures which would satisfy the fixed gas-side temperatures were calculated. The 
fixed heat load and required coolant temperatures, modified by an allowance for pump- 
motor cooling, were used to calculate prime space-radiator area needs (appendix C). 

For the study of simulated solar- or radioisotope- system operation, the 
compressor-inlet temperature was constrained by assigning fixed values for prime 
space- radiator area and coolant mass-flow rate. Turbine- inlet temperature and gas in- 
ventory were also assigned for each calculation. Additional iterations between the gas 
and liquid loops were needed (appendix B). 

RESULTS AND DISCUSSION 

The study results are presented mainly to show trends in conversion efficiency V crw 
to be expected from the system ground-test configuration. Calculations were extended 
beyond design conditions to indicate system potential performance. A promising future 
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change to the compressor-diffuser vane angle was studied for its effect on system oper- 
ating conditions. The results are independent of current hardware limits (see p. 6). 


Parametric System Operation 

Parametric effects on system power output and conversion efficiency were studied 
one at a time. The following parameters were maintained at their design conditions 
(table I) unless they themselves were being studied: turbine- inlet temperature, 
compressor-inlet temperature, working gas, and gas-mixture molecular weight. 
Working-gas inventory, as an independent variable, was studied over a range from 
0. 4 to 1. 6 pounds mass in steps of 0. 4 (0. 18 to 0, 73 kg in steps of 0. 18) for each set of 
assigned parametric values. Results, however, are plotted against gross system output 
power. Gross output includes that part of the system power which would be dissipated 
in the parasitic- load resistors (fig. 1). For initial power system operation, 500 watts 
are maintained in the parasitic load to accommodate variations in system operating 
parameters. For the calculation of space-radiator area needs, an effective radiator 
heat- sink temperature of 450° R (250 K) was assumed. This sink temperature corre- 
sponds to a cylindrical radiator formed around a vehicle whose axis is oriented along 
its path in a low Earth orbit (ref. 7). 

Gas inventory effects . - Figure 2 presents the effects of gas inventory on system 
operating conditions. Increasing gas inventory resulted in increased gas pressure and 
electric power output. The gas- inventory range resulted in from 1.7 to 15. 6 kilowatts 
of output. To maintain the design turbine-inlet temperature, increased heat-source 
input was needed with increased output. 

Figure 2(a) shows the required heat- source, heat- exchanger, ideal thermal input 
Qj, cycle efficiency and conversion efficiency V crw plotted against power output. 

Cycle efficiency was 0. 37 at 2 kilowatts of output power. It increased to about 0. 41 in 
the range from about 10 to 15 kilowatts. The relatively small change in cycle efficiency 
is a result of compensating changes in component performance (explained in more detail 
later in this section). With increasing power output, compressor and turbine efficiency 
and cycle-loss pressure ratio v^/v^ increased, while recuperator effectiveness de- 
creased. The combined positive effect of turbomachinery efficiency and cycle-loss 
pressure ratio exceeded the effect of the reduction in recuperator effectiveness. This 
accounts for the rise in cycle efficiency with increasing output power. 

Conversion efficiency was 0. 17 at 2 kilowatts of output power. It rose to about 0. 30 
in the range from about 10 to 15 kilowatts. Conversion efficiency showed a larger change 
than cycle efficiency because of the added losses from the alternator and system con- 
trols. At the lower power levels, these losses are a larger percentage of the output 
power. Above about 10 kilowatts, conversion efficiency was nearly constant. Conver- 
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Gross electric power output, P S y q, kW 

(a) Thermal input needs for electric power output; cycle and conversion efficiency. 


400 
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(b) Pressure-level changes with gross power 
output 


Figure 2. - Effects of gas inventory on system operation at design conditions. Turbine-inlet temperature, 
2060° R (1144 K); compressor-inlet temperature, 540° R (300 K); working gas, helium-xenon (molecular 
weight, 83.8). 


sion efficiency did not continue to increase because of rapidly increasing alternator 
electromagnetic losses. 

Due to the combined component effects, the required source-heat- exchanger thermal 

input (fig. 2(a)) was not linear with power output. 

Figure 2(b) shows the required compressor-outlet pressure plotted against gross 

power output. For 2 kilowatts of output, a compressor- outlet pressure of 16. 5 psia 
2 

(114 kN/m abs) was needed, while at 15. 7 kilowatts, the required pressure increased 

2 

to 61 psia (421 kN/m abs). Over the range from about 2 to 10 kilowatts, required 
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Turbomachinery Rotating-machinery Effectiveness or 

pressure ratio efficiency pressure ratio 



“ (d-1) Outlet temperatures. 

CD CD 



(c) Changes in rotating-machinery and recu- 
perator performance with gross power output. 



Coolant mass-flow rate, kg/sec 


(d-2) Inlet temperatures. 

(d> Coolant mass-flow rates and temperatures that provide 
the design compressor-inlet temperature. Gas-loop con- 
ditions are constant for each value of gross power. 



Coolant mass-flow rate, Ibm/sec 



Coolant mass-flow rate, kg/sec 


(e) Changes in radiator-area needs with coolant mass-flow 
rate and gross power. Gas-loop conditions are constant 
for each value of gross power output Effective radiator 
heat-sink temperature, 450° R (250 K). 


Figure 2. - Concluded. 
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pressure was almost linear with power level. 

Figure 2(c) shows the predicted rotating-machinery and recuperator operating con- 
ditions plotted against gross system power output. With increasing gross power, the 
required compressor pressure ratio decreased, while the turbine pressure ratio in- 
creased. Cycle-loss pressure ratio, used as a relative measure of total system pres- 
sure losses, increased from 0. 895 at 2 kilowatts to 0. 938 at 15. 6 kilowatts. This trend 
resulted from the heat- exchanger and ducting pressure drops increasing at a lower rate 
than the gas pressure level itself. At 10 kilowatts, the cycle-loss pressure ratio is 
predicted to be 0. 928, rather than the analytical design value of 0. 92 (table I). This was 
a result of conservatism in assigning heat-exchanger and ducting design pressure drops. 

Predicted recuperator heat-transfer effectiveness decreased from about 0. 96 at 
2 kilowatts to about 0. 93 near 16 kilowatts. The recuperator heat-transfer coefficient 
increases with increasing gas mass-flow rate or power output. However, the amount 
of heat to be transferred over the fixed core surface area also increases. The net ef- 
fect is a decrease in effectiveness with increasing power output. At 10 kilowatts, the 
recuperator effectiveness is about 0. 94 (see appendix C) rather than the design value of 
0.95 (table I). 

Rotating-machinery efficiencies are also shown in figure 2(c). Both turbine and 
compressor efficiencies were based on total temperature and pressures. Since the tur- 
bine and compressor are small (ref. 2) and operate over a range of gas densities, vis- 
cous effects on their efficiencies were considered. Reynolds number corrections were 
made on the efficiencies obtained from turbomachinery maps to account for the viscous 
effects (appendix C and ref. 8). The increase in turbine and compressor efficiency 
(fig. 2(c)) with increasing output power were primarily due to the Reynolds number cor- 
rection. The effect was much greater for the turbine than for the compressor. Over 
the power range in figure 2(c), turbine efficiency changed by about 0. 04, while com- 
pressor efficiency changed by about 0. 01. 

At 10 kilowatts of gross system output, turbine efficiency was about 0. 91 and com- 
pressor efficiency was 0.79. The compressor operating point with the design 
compressor-diffuser vane angle was below the peak compressor efficiency of 0. 82. 

Alternator electromagnetic efficiency had a peak value of about 0. 91 at 6 kilowatts 
of gross system output. Inclusion of windage losses reduced the overall alternator ef- 
ficiency. The peak overall alternator efficiency was 0. 86 and occurred at about 7 kilo- 
watts of gross system output. Alternator windage, although reduced in level, was a 
bigger part of the output at the lower power levels. This accounted for the increasing 
difference between electromagnetic and overall alternator efficiency with decreasing 
system output. 

Figure 2(d) shows those combinations of coolant mass-flow rates and waste-heat- 
exchanger coolant temperatures that provide the design compressor-inlet temperature. 
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Matching curves are shown for four gross system output power levels. The power 
levels correspond to four values of constant gas inventory. For fixed gas- loop condi- 
tions, increasing coolant mass-flow rate resulted in increasing heat -transfer effective- 
ness. The coolant-inlet temperature rose asymptotically toward the fixed compressor- 
inlet temperature, while the resulting coolant- outlet temperature dropped because of the 
fixed heat load. As a result of these temperature changes, there was one value of cool- 
ant mass-flow rate for each output power level that resulted in a minimum prime 
radiator area. The heat rejected per unit of radiator surface area varies as the fourth 
power of the radiating temperature of that segment. Hence, there is a trade-off between 
decreasing radiator- inlet temperature and increasing radiator -outlet temperature. 

Figure 2(e) shows the prime radiator area needs plotted against coolant mass-flow 

rate for the four values of gross output power. The locus of minimum prime area 

needs was approximately linear with coolant mass-flow rate. Lines are crossplotted on 

figure 2(d) to show the waste-heat- exchanger coolant temperatures associated with the 

locus of minimum-area mass-flow rates. On the basis of minimum specific prime radi- 

2 

ator area, a ratio A /P „ of about 80 square feet per kilowatt (7.4 m /kW) was 

p, ni sy> ij 

needed at about the 2 -kilowatt level, decreasing to about 46 square feet per kilowatt 
(4. 3 m^/kN) in the range from 10 to 15 kilowatts. Due to simplifying assumptions in the 
calculation of radiator area (see pp. 33 and 55), the actual minimum area needs will be 
about 20 percent larger. 

T urbine- inlet temperature effects. - Turbine- inlet temperature was studied over a 
range from 1660° to 2160° R (923 to 1200 K) in 100° R (56 K) changes. The lowest tem- 
perature was picked to determine conversion efficiency at zirconium-hydride-reactor 
temperature levels. The highest temperature, which exceeds the design value by 100° R 
(56 K), was examined to establish trends. 

Figure 3 shows the effects of varying turbine- inlet temperature on system operation 
at design compressor- inlet temperature. Working-gas inventory was the independent 
variable. Figure 3(a) shows conversion efficiency and minimum prime radiator area 
needs plotted against power output. (In contrast to fig. 2(d), only the minimum areas 
are shown in fig. 3(a).) Lines of constant source- heat- exchanger ideal thermal input are 
shown, and lines of constant conversion efficiency are crossplotted on the radiator area 
curves. Figure 3(b) shows compressor-outlet pressure plotted against gross power 
output. Lines of constant gas inventory and conversion efficiency are crossplotted on 
the curves. 

Trends in system operating conditions at each value of turbine- inlet temperature are 
similar to those at the design value. Maintaining constant gross system power output 
with decreasing turbine- inlet temperature requires increased source thermal input 
(fig. 3(a)) and increased gas inventory and compressor-outlet pressure (fig. 3(b)). The 
resulting conversion efficiency decreases, while the required minimum prime radiator 
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Source-heat-exchanger 
ideal thermal input, 



(a) Changes in efficiency and radiator area needs with gross 
power output and turbine-inlet temperature. Effective 
radiator heat-sink temperature, 450° R (250 K). 



0 4 8 12 16 20 

Gross electric power output, P sy G , kW 


(b) Changes in compressor-outlet pressure with gross 
power output and turbine-inlet temperature. 

Figure 3. - Effects of turbine-inlet temperature on system 
operation at design compressor-inlet temperature. Work- 
ing gas, helium-xenon (molecular weight, 83.8). 
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areas increase (fig. 3(a)). Maintaining constant source thermal input (fig. 3(a)) or gas 
inventory (fig. 3(b)) with decreasing turbine- inlet temperature results in decreasing sys- 
tem power output. 

System operation at 2160° R (1200 K), if it is possible, would gain about 0. 015 in 
conversion efficiency over that at the design inlet temperature. At a turbine-inlet tem- 
perature of 1660° R (923 K), the highest studied value of gas inventory (fig. 3(b), 

1. 6 lbm or 0. 73 kg) resulted in a gross system power output of 7. 5 kilowatts. Conver- 
sion efficiency at 1600° R increased from 0. 11 at 2 kilowatts to about 0. 21 at 7. 5 kilo- 
watts. 

Compressor- inlet temperature effects . - Compressor- inlet temperature was studied 
over a range from 520° to 600° R (289 to 333 K) in 20° R (11 K) changes. The upper 
value was chosen to be typical of reactor-type system operation, while the lower value 
was used to establish potential conversion efficiency gains over those at the design inlet 
temperature. 

Figure 4 shows the effects of varying compressor- inlet temperature on system op- 
eration at design turbine-inlet temperature. The same format as in figure 3 was re- 
peated in figure 4. 

Trends in system operating conditions at each value of compressor -inlet tempera- 
ture are similar to those at the design value. Maintaining constant power output with 
increasing compressor-inlet temperature requires increased source thermal input 
(fig. 4(a)) and increased gas inventory and compressor-outlet pressure (fig. 4(b)). The 
resulting conversion efficiency and required minimum prime radiator area both decrease 
(fig. 4(a)). Maintaining source thermal input (fig. 4(a)) or gas inventory (fig. 4(b)) con- 
stant with increasing compressor-inlet temperature results in decreasing system power 
output. 

System operation at 520° R (289 K) would gain about 0. 01 in conversion efficiency 
over that at the design inlet temperature; however, the required radiator area would in- 
crease. At a compressor- inlet temperature of 600° R (333 K), the highest studied value 
of gas inventory (fig. 4(b), 1. 6 lbm or 0. 73 kg) resulted in a gross system power output 
of 12. 3 kilowatts. Conversion efficiency at 600° R increased from about 0. 15 at 2 kilo- 
watts to about 0. 26 in the 9- to 12-kilowatt range. 

Operation with krypton. - Krypton has a molecular weight of 83. 8, and its viscosity 
is about the same as the design gas mixture. Therefore, the fluid-dynamic losses with 
either krypton or the design mixture are about the same. Krypton has a thermal con- 
ductivity lower than that of the mixture, which reduces heat transfer. Initial system 
tests will use krypton as the working gas since it is much less expensive. 

Comparative system operation with krypton and the design mixture of helium and 
xenon is shown in figure 5 for the design turbine- and compressor-inlet temperatures. 
Conversion efficiency and minimum required prime radiator area are plotted against 


15 










Gross electric power output, P G , kW 

Figure 5. - Comparative system operation with 
krypton and design gas mixture. Turbine-inlet 
temperature, 2060° R (1144 K); compressor- 
inlet temperature, 540° R (300 K); effective 
radiator heat -sink temperature, 450° R (250 K). 

gross system output power. Four lines of constant source ideal thermal input are cross- 
plotted on the results. 

Due to the lower thermal conductivity of krypton, conversion efficiencies were 
lower, and minimum required prime radiator areas were larger, than those with the use 
of the gas mixture. More heat-source thermal input was needed to obtain the same out- 
put power obtained with the design mixture. The same compressor-outlet pressures 
were required for both the mixture and krypton at any value of gross output power. With 
the use of krypton, conversion efficiency increased from about 0. 16 at 2 kilowatts of out- 
put to about 0. 27 in the range from 10 to 16 kilowatts. Compared with the design gas 
mixture, conversion efficiency was lower by 0. 01 at the 2-kilowatt level and by about 
0. 03 in the 10- to 16-kilowatt range. 

Gas-mixture relative helium mass effects. - For the design molecular weight of 
83. 8, the helium mass is 1.78 percent of the total gas mass. The helium volume is 
37. 3 percent of the total gas volume. To determine the sensitivity of system conversion 
efficiency due to errors in gas mixing, the effects of ±20 percent changes in the relative 
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helium mass were studied at the design turbine- and compressor-inlet temperatures. 

At all system gross output power levels in the range from 2 to about 16 kilowatts, 
chatiges in conversion efficiency were small with the ±20 percent change in relative 
helium mass. The largest changes occurred at about the 8-kilowatt output power level. 
There, the ±20 percent mass change caused a ±0. 004 change in conversion efficiency. 
With increased helium mass content, the compressor-map operating point, although at 
a lower corrected speed, shifted more toward peak efficiency. Increasing helium mass 
content increased compressor efficiency and heat-transfer effectiveness, but decreased 
turbine efficiency. 

Further increases in relative helium mass above 20 percent would result in addi- 
tional small gains in conversion efficiency. At some point the trend would reverse. 
However, determining the relative helium mass increase needed to maximize conversion 
efficiency was beyond the purpose of this study. For errors in gas mixing as large as 
±20 percent of the relative helium mass, there would be little effect on system conver- 
sion efficiency. 


Simulated Solar Operation 

Calculations were made for simulated solar-power- system operation. Correspond- 
ing to the source configuration in reference 2, a solar-heat receiver designed for 40. 5 
kilowatts of average thermal input to the design gas mixture at a turbine- inlet tempera- 
ture of 1960° R (1089 K) was assumed. Operation was assumed with a simulated prime 
radiator area of 500 square feet (46. 5 m^). An effective radiator heat- sink tempera- 
ture of 400° R (222 K) was used. This sink temperature corresponds to a cylindrical 
radiating surface that is oriented parallel to the direction of the Sun and in a low Earth 
orbit (ref. 9). Since the solar collector must be oriented with respect to the Sun, the 
lower sink temperature for this case is reasonable. 

The study was made to predict changes in conversion efficiency with changes in gas 
inventory for this type of ground- simulated, fixed- system operation. Turbine- inlet tem- 
perature was kept at 1960° R (1089 K) independent of gas inventory. In space, with a 
fixed solar-heat receiver, turbine-inlet temperature would change somewhat with gas 
inventory and orbital position. 

For any fixed radiator area and design level of source thermal input, there is a 
particular value of coolant mass-flow rate that maximizes conversion efficiency and 
power output. This optimum coolant mass-flow rate also results in the lowest value of 
compressor-inlet temperature. For the calculations presented herein, the optimized 
coolant mass-flow rate was 0. 25 pound mass per second (0. 11 kg/sec). 
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Gross electric power output, P S y q, kW 


(a) Thermal input needs and compressor- (b) Efficiency and pressure changes 

inlet temperature changes with gross with gross power output, 

power output. 


Figure 6. - Simulated solar operation. Turbine-inlet temperature, 1960° R (1089 K); working gas, 
helium-xenon (molecular weight, 83.81; effective radiator heat-sink temperature, 400° R <222 K); 
prime radiator area, 500 square feet (46. 5 m^l; coolant mass-flow rate, 0. 25 pound mass per sec- 
ond (0. 11 kg/sec); design average thermal input, 40. 5 kilowatts. 


Figure 6 shows the changes in simulated solar- power- system operation due to 
changes in working-gas inventory. The optimized coolant mass-flow rate was used. 
Inventory was varied from 0. 8 to 1.6 pounds mass in steps of 0. 2 (0. 36 to 0. 73 kg in 
steps of 0. 09). However, the results are plotted against gross system power output. 
Figure 6(a) shows the variations in compressor-inlet temperature and required source 
ideal thermal input, while figure 6(b) shows conversion efficiency and required 
compressor-outlet pressure. 

The design average thermal input of 40. 5 kilowatts resulted in 11.7 kilowatts of 
gross system power output. At 500 square feet of prime radiator area, 11.7 kilowatts 
of output was the maximum. The corresponding operating conditions were 
(1) compressor-inlet temperature, 526° R (292 K), (2) compressor-outlet pressure, 

48.7 psia (336 kN/m 2 abs), and (3) conversion efficiency, 0.29. 

For the simulated fixed-system operation, increasing gas inventory resulted in in- 
creasing power output. Compressor-inlet temperature, required source ideal thermal 
input, and compressor-outlet pressure all increased with increasing gas inventory. 
System conversion efficiency showed small changes with output power. At 7. 1 kilowatts, 
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conversion efficiency was 0. 292, while at 13. 1 kilowatts it was 0. 278. Conversion ef- 
ficiency reached a peak value of 0. 298 at 10 kilowatts. 

For any fixed coolant mass-flow rate and radiator area, there are optimum values 
of source thermal input power and gas inventory that correspond with a peak in con- 
version efficiency. The value of input power at which peak efficiency occurs is pri- 
marily a function of the radiator area. The larger the radiator area, the higher the 
value of thermal input at peak efficiency. In figure 6(b), the difference between the con- 
version efficiency at 11.7 kilowatts and the peak efficiency at 10 kilowatts is small, 
about 0. 01. A reoptimization in both coolant mass-flow rate and radiator area would 
result in only a small efficiency gain for the design average thermal input power level. 

This calculation of simulated solar-power-system operation has shown that a space 
version with a design thermal power level and a limited amount of available radiator 
area should be optimized with coolant mass-flow rate to provide maximum power output. 
If radiator area is not limited, it could be chosen to prefer either design or off-design 
thermal input power operation; that is, a user could choose a large radiator area to 
maximize power output at design thermal input power, or he could pick a smaller area 
to keep a relatively high conversion efficiency over a range of lower-than-design input 
power. 


Simulated Radioisotope Operation 

Calculations were made for simulated radioisotope- power- system operation. Cor- 
responding to the heat source in reference 2, a radioisotope source rated at 25 kilowatts 
thermal was studied. All of this thermal energy was assumed to radiate to the heat ex- 
changer and be absorbed without loss by the system working gas. In a space version, both 
the radioisotope capsule array and adjacent source heat exchanger would be sized to 
provide the 25-kilowatt-thermal input to the gas at a turbine-inlet temperature of 
2060° R (1144 K). Operation was assumed with a simulated prime radiator area of 
350 square feet (32. 5 m). An effective radiator heat-sink temperature of 450 R 
(250 K) was used. The coolant mass-flow rate, optimized for the assumed radiator area, 
the heat-source power level, and the design turbine- inlet temperature, was 0. 155 pound 
mass per second (0. 704 kg/sec). This value was held constant in the study. 

The study was made to predict changes in conversion efficiency with changes in both 
working-gas inventory and radioisotope thermal power rating for the fixed-power-system 
operation. Since computed results give the required source ideal thermal input for an 
assigned turbine- inlet temperature, crossplotting was needed to obtain the results for 
constant values of source ideal thermal input. 
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(a) Inlet temperature changes with gross (b) Efficiency and pressure changes 

power output and source ideal thermal with gross power and source ideal 

input. thermal input. 


Figure 7. - Simulated radioisotope operation. Working gas, helium-xenon (molecular weight, 33.8); 
effective radiator heat-sink temperature, 450° R (250 K); prime radiator area, 350 square feet 
(32. 5 m 2 ); coolant mass-flow rate, 0. 155 pound mass per second (0. 704 kg/sec). 


Figure 7 shows the changes in simulated, radioisotope- power- system operation due 
to changes in both working-gas inventory and source ideal thermal input power. The 
results for the assigned radiator conditions and the design gas mixture are plotted 
against gross system output power. Solid curves are used to present the results for 
three constant values of source-heat-exchanger ideal thermal input: 20, 25, and 30 kilo- 
watts. Figure 7(a) shows the values of turbine- and compressor- inlet temperature. 
Figure 7(b) shows the values of conversion efficiency and compressor-outlet pressure. 

Figure 7 shows that for a constant value of source thermal input power, an initially 
high value of gas inventory would result in the following relative operating conditions: 
high compressor-outlet pressure, low turbine- inlet temperature, high compressor-inlet 
temperature, and low conversion efficiency. Bleeding gas from the system would then 
improve operation. 
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At design turbine- inlet temperature with an ideal 25-kilowatt, source thermal input, 

2 

the compressor-outlet pressure was 31. 3 psia (216 kN/m abs). The corresponding 
compressor -inlet temperature was 530° R (284 K), gross power output was about 7 kilo- 
watts, and conversion efficiency was about 0. 28. 

Figure 7(b) shows that, at the design turbine- inlet temperature conditions, a maxi- 
mum conversion efficiency occurred for about 25 kilowatts of ideal thermal input. As 
was the case for simulated solar-power-system operation, with fixed values of radiator 
area and coolant mass-flow rate, there were particular values of source input and gas 
inventory that maximized power output for operation at the design turbine-inlet tempera- 
ture. 

The turbine- inlet-temperature curves (fig. 7(a)) for 25 and 30 kilowatts of thermal 
input intersect at about 4. 6 kilowatts of output. For constant, but lower-than-design, 
turbine- inlet temperatures, the same fixed radiator conditions also result in unique 
maximum power output conditions. At reduced turbine- inlet temperatures, both the out- 
put power level and the value of source thermal input for maximum conversion efficiency 
operation are reduced. At a turbine-inlet temperature of 1740° R (966 K), the change in 
conditions for maximum power output was enough that increasing source input above 
25 kilowatts did not result in more output power. 

The constant-design-turbine-inlet-temperature curves (long-dashed curves) of fig- 
ure 7 also show a way in which power output could be controlled during decay of a short- 
lived radioisotope. If a source was rated for 30 kilowatts input at the "start of 
mission, " the system could initially produce about 8. 5 kilowatts of gross power at the 
rated turbine-inlet temperature. Conversion efficiency, about 0. 28, would be only 
slightly below its peak value. As the source decays, and if no change is made to the 
system operation, its power output would drop as indicated by the constant-gas- inventory 
line (short-dashed curve) in figure 7(b). After decay to 25 kilowatts (a 16.7 percent 
decrease), the gross output power would drop to about 5. 5 kilowatts (or a 36 percent 
decrease). However, if the working gas is bled off (reducing system pressure level) to 
maintain the design turbine-inlet temperature during source decay from 30 to 25 kilo- 
watts, gross power output would drop to 7. 1 kilowatts, or only 16 percent. From the 
upper curves in figure 7(b), such operation would maintain a relatively high conversion 
efficiency over a 33 percent decay from the initial 30 kilowatts. 


Reset Compressor-Diffuser Vanes 

Up to this section all results were based on the present hardware intended for use 
in the system ground tests. The design compressor-diffuser vane angle was chosen to 
provide margin from compressor surge as protection in case of mismatched system op- 
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erating conditions. This choice resulted in a compressor-corrected mass-flow rate, at 
system design conditions, above that corresponding to peak compressor efficiency. Re- 
sults of unpublished tests on this compressor at Lewis, over a range in vane angles, 
have shown that a 3° increase would result in a shift to peak compressor efficiency for 
design system operation. 

This analysis has shown no mismatch in system operating conditions. Hence, if this 
is confirmed in the system ground tests, the use of a modified, +3° compressor -diffuser 
would be desirable. 

Figure 8 shows the potential system operational gains with the use of a +3° 
compressor-diffuser. Results are shown for both the design and modified compressor- 
diffuser at the design turbine- and compressor- inlet temperatures plotted against gross 
system output power. The design gas mixture was used and the range in gas inventory 
from 0. 4 to 1. 6 pounds mass (0. 18 to 0.73 kg) was studied. Figure 8(a) shows conver- 
sion efficiency and minimum required prime radiator area, while figure 8(b) shows the 
required compressor-outlet pressures. Four lines of constant source ideal thermal 
input are crossplotted on the results in figure 8(a). The four lines of constant gas in- 
ventory are shown in figure 8(b). 

With the modified compressor-diffuser, conversion efficiency increased from about 
0. 18 at 2 kilowatts of output to about 0. 32 in the 10- to 15-kilowatt power range. Over 
the 2- to 15-kilowatt output range, the gain in conversion efficiency with the modified 
compressor was between 0. 01 and 0. 02. The gains in conversion efficiency were pri- 
marily due to increased compressor efficiency. At the 10-kilowatt level with the +3 U dif- 
fuser, compressor efficiency was at its peak value of about 0. 83. (The modified diffuser 
increased the peak compressor efficiency by about 0. 01 over that of the design diffuser, 
see p. 12). 

More power output for the same thermal power input was obtained with the use of 

the modified diffuser. Over the power output range, use of the +3° diffuser reduced the 

minimum required specific prime radiator area by between 4 and 5 square feet per kilo- 
o 

watt (0. 4 and 0. 5 m /kW). Figure 8(b) shows that, for the same output power, the use 
of the modified diffuser resulted in lower required compressor-outlet pressure and gas 
inventory. 


CONCLUDING REMARKS 

This analysis approximated the power system configuration being ground tested. 
Because some real- system effects were neglected (see p. 7), the computed levels of 
system conversion efficiency must be considered ideal. However, the trends of the re- 
sults should be correct. 
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The neglected effects may lower the level of system conversion efficiency shown in 
the analysis by about 20 percent. Therefore, before experimental results are available, 
system conversion efficiency at design conditions should be considered to be about 0. 14 
at 2 kilowatts of output, increasing to about 0. 24 in the 10- to 15-kilowatt range. All 
other trends shown in the analytical results should be referenced to this reduced effi- 
ciency level. 


SUMMARY OF RESULTS 

A steady- state digital- computer study of a Brayton power system was made to ex- 
amine parametric and simulated-fixed- system operations to be conducted in experimen- 
tal system ground tests. The analysis shows trends in conversion efficiency and is in- 
dependent of present system operational limits. However, the present limitations of 
the existing hardware are discussed. 

Unless noted otherwise, the results are summarized for the design working gas, 
a mixture of helium and xenon at a molecular weight of 83. 8; the design turbine-inlet 
temperature, 2060° R (1144 K); and the design compressor -inlet temperature, 540° R 
(300 K). Major study results were as follows: 

1. Estimated levels of system conversion efficiency, including an allowance for 
losses neglected in the analysis, increased from about 0. 14 at 2 kilowatts of output to 
about 0. 24 in the 10- to 15-kilowatt range. The ideal, predicted levels were 0. 17 at 
2 kilowatts and 0. 30 in the 10- to 15-kilowatt range. 

2. Future power system gains were predicted with the use of a compressor-diffuser 
modified by a 3° increase in vane angle. With the modified diffuser, the analysis 
showed potential conversion efficiency gains of 0. 01 at the 2-kilowatt output power level, 
and of 0. 02 in the 10- to 15-kilowatt range. 

3. Calculations for both simulated solar- and radioisotope-power- system operation 
showed that with fixed values of space radiator area and coolant mass-flow rate, there 
were particular values of source thermal input power and gas inventory that maxi- 
mized gross system power output. 

4. Simulated radioisotope-power- system calculations showed that bleeding working 
gas would be an effective means to maintain high system electric power output during 
the decay of a short-lived radioisotope. 

5. Calculations showed that, with the use of krypton as the system working gas, 
conversion efficiency would decrease by about 0. 01 at 2 kilowatts of output and by about 
0. 03 in the 10- to 15-kilowatt range. 
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6. A ±20 percent change in the relative helium mass of the working-gas mixture 
caused no more than a ±0. 004 change in conversion efficiency within the 2- to 15-kilowatt 
output power range. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 25, 1969, 

120-27. 
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APPENDIX A 


SYMBOLS 


A 

A,. 


w 


p 

D t 

E 

f 

f',g',h' 


f 


cw 


G 

g 

h 

i 

J 

3 

K 

K 

eq 

k 

Lw 


9 9 

heat-transfer surface area, ft"'; m"" 

2 2 

ideal prime space radiator area (see p. 33 and eq. (C103)), ft ; m 
waste-heat- exchanger capacity-rate ratio, mc p y/(mCp^ 

constant defined in eq. (C9), (°R)(lbm)/(ft)(lbf)(rpm); (K)(kg)/(N)(m) 

constant defined in eq. (C21), dimensionless 

specific heat at constant pressure, Btu/(lbm)(°R); J/(kg)(K) 

tip diameter, ft; m 

heat-transfer effectiveness 

friction coefficient 

functions of ( ) 


compressor work factor, 


(60)%Jc p (T 5 - T 4 ) 

" 2n X c 


dimensionless; 


gc n ( T 5 - TO 4 
p “ ■ 2 * 4 X 10* 

2 xt 2t,2 

7T N D t, c 

2 2 

mass-flow rate per unit area, lbm/(ft )(sec); kg/(m )(sec) 
acceleration due to gravity, 32.2 (ft)(lbm)/(sec^)(lbf); 9.81 (m)(kg)/(sec^)(N) 
local heat-transfer coefficient, Btu/(sec)(ft^)(°R); J/(sec)(m^)(K) 
a number index 

mechanical equivalent of heat, 778 (ft) (lbf)/Btu; 1 (m)(N)/J 

Colburn factor for heat transfer, (h/Gc p )(Pr)^ // ^, dimensionless 

molecular weight correction factor, 83. 8/M 

w 

equivalent duct length-to-diameter ratio (eq. (C94)), dimensionless 
thermal conductivity, Btu/(sec)(ft)(°R); J/(sec)(m)(K) 
bearing friction loss, kW 
alternator windage loss, kW 
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M, 

m 

N 


w 


N 


s, C 


tu 


N, 

n 

n 

P 

P 

P 

e 

Pr 


sy, G 


Qj 

Re 




T'‘C 


T 

T £ 

U 

V 

v 


molecular weight 

mass-flow rate, Ibm/sec; kg/sec 

shaft rotational speed, rpm; rad/sec 


compressor specific speed, 



1800 


eJc T frfr" 1 )/’' 
g c p 4\ C 



dimensionless; 



S c p T 4 


{y-l)/y 



number of thermal units, UA 



dimensionless 


empirical exponent 
number of passes 
power, kW 

internal system electric power needs, kW 
Prandtl number, pCp/k, dimensionless 
gross system electric power output (fig. 1), kW 
total gas pressure, psia; kN/m abs 
thermal power, kW 

ideal thermal power input to gas in the heat-source heat exchanger, kW 

Reynolds number, pvD|./ p, dimensionless 

compressor total-pressure ratio, P5/P4 

turbine total- pressure ratio, P4/P2 

cycle-loss pressure ratio 

total temperature, °R; K 

effective radiator heat- sink temperature, °R; K 

overall heat-transfer coefficient, Btu/(sec)(ft^)(°R); J/(sec)(m^)(K) 

3 3 

volume, ft ; m 
velocity, ft/sec; m/sec 
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W gas inventory, lbm; kg 

r mathematical function defined in eq. (C78) 

y ratio of specific heats for a monatomic ideal gas, 5/3 

A difference operator 

6 ratio of compressor-inlet pressure to U. S. standard pressure, p 4 /14. 69 psi; 

P 4 /101. 325 kN/m 2 

e radiator surface hemispherical emittance, 0 . 90 

£ mathematical function defined in eq. (C80) 

77 efficiency 

^cnv conversion efficiency, P gy 

V cy cyc le efficiency, P g ^ G /Q. 

6 ratio of compressor- inlet temperature to U. S. standard temperature, 

T 4 /519 °R; T 4 /288 K 

p viscosity, lbm/(ft)(sec); kg/ (m) (sec) 

p density, lbm/ft 3 ; kg/m 3 

a Stefan-Boltzmann constant, 0. 481xl0“ 42 Btu/(sec)(ft 2 )(°R 4 ); 5.73X10 ” 8 

J/(sec) (m 2 ) (K 4 ) 

r turbine torque, ft/lbf; m/N 

cp waste-heat- exchanger performance factor defined in eq. (C76) 

Subscripts: 
a alternator 

C compressor 

c cold side of recuperator 

cs compressor surge 

em electromagnetic 

f fin 

G gross 

g gas 

h hot side of recuperator 

i inlet 
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I 


liquid 
m minimum 

n net 

0 outlet 

P coolant pump 

p per pass 

R radiator 

r recuperator 

s surface 

sh shaft 

so source heat exchanger 

T turbine 

w waste heat exchanger 

1 turbine- inlet station 

2 turbine- outlet station 

3 waste-heat- exchanger, gas- side, inlet station 

4 compressor- inlet station 

5 compressor- outlet station 

6 source-heat- exchanger inlet station 

Superscripts: 

* design or reference condition 

average value 
" new value 
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APPENDIX B 


DETAILS OF COMPUTER PROCEDURES 
Fluid Properties 


Values of gas viscosity and thermal conductivity were tabulated as functions of tem- 
perature and entered into the computer program as input data. These properties were 
obtained from Vanco's analysis (ref. 10). Pressure effects on the properties were small 
for the range of interest and were neglected. In the computer program, all tabular in- 
terpolations used Newton's formulas (ref. 11, pp. 61-66). 

Coolant specific-heat values used in the liquid- loop calculations were obtained from 
testing three sample batches of the coolant over the temperature range from 77° to 
300° F (298 to 422 K). The following linear fit was found and used in the program: 

c j = 3. 10 T Xl0~ 4 + 0.239 

P, 6 l 

c = 2. 34 T + 1000 
Pj ^ v 

(Where two equations are given, the first is in U. S. customary units and the second in 
SI units. ) 


(Bl) 


Flow-Logic Diagrams 

Flow-logic diagrams for the computer calculations are shown in figures 9 and 10. 
Figure 9 shows the logic paths for gas-loop closure which is common to both types of 
calculations in figure 10. Figure 10(a) shows the logic paths for the calculation of sys- 
tem parametric operation. And figure 10(b) shows the logic paths for simulated-fixed- 
system -operation calculations. 

Gas-loop closure . - For either calculation, this part of the computer program is 
entered with assigned values for W_, M , N, T*, and T 4 . Values are assumed for 

gw 1 *x 

m and pj. During the first pass (n p = 1) through this calculation, all ratio terms for 
heat- exchanger pressure drops and recuperator effectiveness, other than (m/m*) n , are 
set equal to 1. The logic requires more than one pass for each new entry. After the 
first pass, turbine- inlet temperature is lowered because of the seal loss model, and p^ 
is corrected by -(pj’ - Pj)/2. 
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Figure 9. - Flow-logic diagram for gas-loop closure. 



























(a) Calculation of system parametric operation. 



(b) Calculation of simulated fixed-system operation. 


Figure 10. - Flow-logic diagrams for calculation of system operation. 


The first iteration to be satisfied is that shown for turbine- inlet pressure Pj. The 


second iteration, in series, is on working-gas inventory, 
compared to the assigned value W in the problem input. 

D 


Calculated inventory W^' is 
If gas inventory is not satis- 


fied, both m and p^ are corrected and the process repeats. The changes are made 


based on linear variations of m and p. with W . Both m and p. are corrected for 

1 g 1 

_(W” - Wg)/2 times the respective linear ratio. 

When both p^ and are satisfied, this part of the calculation is finished 
and the calculations shown in figure 10 continue. 
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Parametric system operation . - input values (fig. 10(a)) for this calculation include 

W , M w , N, Tj, and T^. After the gas-loop calculations are satisfied, the waste- 

heat- exchanger capacity-rate ratio C is varied over a range of values from 0. 5 to 

w 

1. 0 in steps of 0. 05. The liquid mass-flow rate is calculated from the liquid capacity 
rate (me \ by evaluating c , from equations (Bl) at T , and dividing. After the 

\ P/£ Pj t VV , i 

range in C is complete, the remaining calculations for gross electric power output 
w 

P ~ are made, and the results are printed out. 
sy, G 

Simulated fixed- system operation. - This type of calculation (fig. 10(b)) is a varia- 
tion on the previous methods (fig. 10(a)). Here, A p and m^ are assigned in place of 
T 4 in the input values. Initially, an arbitrary value of T^ is assumed, and the gas- 
loop closure calculations are made. The waste-heat- exchanger calculations are made 

with an iteration on the temperature at which c 7 is evaluated. 

p, t o 

The assigned prime radiator area is divided into 10- square- foot (0. 93-m ) sections. 
Each section is assumed to be isothermal and to reject heat radiantly at its liquid- inlet 
temperature. Hence, the calculation is approximate. It neglects the temperature drops 
across the section and between the coolant and the radiating surface. The inlet tem- 
perature to the first radiator section is Tg Then, in sequence, the outlet tempera- 
ture of each section is calculated from a heat balance, and is the inlet temperature for 
the next radiator section. The process ends when the outlet temperature of the last sec- 
tion T p q is obtained. The coolant pump outlet temperature results from (see 
eq. (Cl 02)). 


(B2) 


If the liquid- loop temperatures T p Q and T^ ^ are not satisfied, the wrong T^ 
was used. Then, T^ is corrected by (T^ ^ - T p q )/2 and the calculation repeats. 
After satisfying both the gas- and liquid- loop temperatures, the remaining calculations 
for gross output power are made, and the results are printed out. 


P, o ~ T R, o + 
P, o = T R, o + 
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APPENDIX C 


COMPONENT MODELS 
Rotating Machinery 

Gross shaft power P gh G was defined as the difference between turbine work and 
compressor work. The net shaft power was then 


P sh, n P sh, G ~ H) ~ Sv 
and the gross alternator power was 


P a, G ^em P sh, n 


(C2) 


Turbine . - The turbine performance maps were obtained from Nusbaum and 
Kofskey's results (ref. 3). Tables entered into the computer program were of the form 


(C3) 



and 



(C4) 


For each turbine calculation, the tables were entered with known values of N/^/kTj and 
m|/^T^/pj. A double interpolation, using Newton's formulas, calculated r T for these 
values. Then another double interpolation solved for t/p^. 

Monatomic perfect-gas relations were assumed, thus 


r = 


(C5) 
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and 


y / L 54 s\ 
Cp “ y - 1 \JM w j 

c = _^/M22\ 




> 


Turbine efficiency resulted from 


Tjrp ~ 


(—) 

ImKTj 


1 - 


r -(r-i)/r 

T 


where 


tN 

mKT l m|/KT[ 

Pi 



(C6) 


(C7) 


(C8) 


and 


C 1 =— — = 2.27x10 
1 30Jc* 

Ir 


C* = — = 4. 02x10" 3 

r * 

P 




> 


(C9) 


with c* given by equations (C6) evaluated for M = 83. 8. A Reynolds number correc- 

p w 

tion for efficiency was assumed, using the general form (see Shepherd, ref. 8) 


A (t?) = (1 - v*) 




(CIO) 
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Nusbaum and Kofskey's results were used to adapt equation (CIO) for use with the turbine 
maps. Turbine Reynolds number was defined as 


Re” 


2m 



T 


This number is related to the conventional usage by 

Re” = - Re T 

7 T 


(Cll) 


(Cl 2) 


The resulting correction was 


7. 62x1 0 4 

Re” 

^ 1 

such that 

7 y” = ?7 T + A(tj t ) (Cl 4) 

Compressor . - The compressor experimental performance results have not been 
published. Figure 11, however, is a sketch of the compressor map formats that were 
used. A transformation was made before entering tabular values into the computer pro- 
gram for ease of interpolation with Newton's formulas. All compressor input data were 


vO.25 


(Cl 3) 


A(7] t ) = 0. 112 




Compressor corrected mass-flow rate, mVe/6 
Figure 11. - Experimental compressor maps. 
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referenced to the estimated compressor surge lines, shown dashed in figure 11. The 
computer tables were of the form 


where 


and 



(Cl 5) 


(Cl 6) 


(C17) 


(C18) 


(Cl 9) 


For each compressor calculation, the tables were entered with known values of N/VK0 
and m\/K0/5. A single interpolation and equation (C18) resulted in A (m\Zfl/6). Two 
double interpolations gave values for r^ and f . 
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Compressor efficiency resulted from 



where 


(C20) 


»Xc| 


C 2 S 


N 

Ve, 


\*2 




519(25)gJc* 


= 0. 556 


>- 


D. 


t,c 


_N_ 


\*2 


288gc* 

r 


y 


(C21) 


with c* given by equations (C6) evaluated for M = 83. 8. The same Reynolds number 
3? w 

correction form (eq. (CIO)) was used for the compressor. And, similar to the defini- 
tion for turbine Reynolds number, compressor Reynolds number was defined as 


Re C 


2m 

^4 D t, C 


(C22) 


Based on the unpublished compressor performance results, the Reynolds number cor- 
rection was 


5 \0. 06 | 
14x10 5 \ 

Re C / . 

such that 

Vq = V c + &(v c ) 


A(77 c ) = 0. 207 


1 - 


(C23) 


(C24) 
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Alternator . - The alternator electromagnetic efficiency was obtained from Repas and 
Edkin's results (ref. 4). They show this efficiency plotted against gross alternator power 
P „ for power factors of 0. 75 and i. 0. Although a power factor higher than 0. 75 may 

VJ 

be expected in most uses, the 0. 75 results were used here for conservatism. Electro- 
magnetic efficiency was plotted, with use of equation (C2), against net shaft power P gh n 
and extrapolated from 15 to 24 kilowatts. This extension was needed to examine opera- 
tion near 16 kilowatts of gross power P p . The efficiency values were tabulated and 
entered into the computer program in the form 


^a, em f ^ P sh, n^ 


(C25) 


Seal losses . - Based on the analytical design of the rotating machinery, 2 percent of 
the compressor-outlet flow must be bled off to supply pressure to the gas bearings. The 
bleed flow is ducted from the compressor scroll into the bearing and shaft cavities. 
One-half of this flow is estimated to leak through the turbine- end, labyrinth shaft seal. 
This flow reenters the main flow just upstream of the turbine impeller. The other half 
of the bleed flow leaks through the compressor- end seal. It reenters the main flow just 
downstream of the compressor impeller and ahead of the diffuser vanes. The bleed flow 
is at or near the compressor-outlet total pressure. There is a sufficient pressure dif- 
ference for the flow to reenter the main stream at the compressor-impeller exit, since 
there the flow is controlled by the local static pressure. 

The leakage flow into the turbine causes a mixing loss. The leakage flow circulating 
around the compressor vanes causes an increased pumping loss. In the computer pro- 
gram, a difference in gas mass-flow rate was calculated as 


m c 


0.98 

0.99 


m. 


(C26) 


where m c is the gas mass-flow rate from the compressor-scroll outlet to the turbine- 
scroll inlet. This mass-flow-rate difference was accounted for in the recuperator cal- 
culations (eq. (C36)). The bleed flow was also assumed to reduce the effective turbine- 
inlet temperature by 


A(Tj) = 


0. 01 (Tj - T 5 ) 
0.99 


(C27) 
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such that 


T'j = Tj - AfTj) 


(C 28 ) 


The author feels that this penalty (eq. (C27)) is conservative and should compensate for 
both bleed-flow penalties. 

Windage loss. - Alternator windage was assumed to vary as (see ref. 8), 




(C29) 


Then, for an ideal gas, 


P 


oc 


M p 
T 


and, for fixed geometry, 


Re cc 


m 


M 


Hence, assuming constant- temperature operation, or T* = T, 


L w _ /m \ n /g*\ n / M wVpY» \ 5 


and for choked flow, m cc p- therefore, 



(C30) 


(C31) 


When equation (C31) and design predictions for windage losses were used, n = -0. 40 was 
found to give an acceptable fit. The design predictions for windage losses were in 
agreement with later test results (ref. 4). The reference or starred conditions in equa- 
tion (C31) were evaluated at the 10-kilowatt gross output conditions and for operation at 
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860° R (480 K). Hence, the expression was 







> (C32) 


Bearing friction. - Both journal- and thrust-bearing friction losses were assumed 
to vary as (ref. 12) 


^ ,/jAM 2 
Lg \»* n **) 


(C33) 


Hence, this approximation assumed no change in total bearing loss with pressure or 
power level. The author believes this to be a reasonable first-order approximation. 
Based on design estimates, Lg = 0. 27 kilowatt. And j±* was evaluated at 860° R 
(480 K) for the design gas mixture. Therefore, 


1^ = 0. 27/ IL 


\2. 49x10 


N 

5/136 000 j 


Lh = 0.27 IL 


(3.71x10 


N 

-5/13770 j 


A 


> (C34) 


Since N was constant for the results presented herein, the bearing friction loss 
changed only with a change from the design gas mixture. 


Heat -Transfer Models 

Approximate relations were developed for the recuperator and waste heat exchanger. 
Empirical exponents were found by fitting these models to design predictions. The 
models were then modified to reflect changes based on acceptance test data. For the 
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results presented in this report, it was assumed that the source heat exchanger could 
be run hot enough to attain all conditions studied. 

Recuperator. - Recuperator heat-transfer effectiveness is a function of the gas 
mass-flow rate and the heat to be transferred. For the purpose of constructing a model, 
the difference in gas mass-flow rates (eq. (C26)) was neglected. After E f was deter- 
mined, the recuperator -outlet temperatures, including the flow difference, were calcu- 
lated from 


Tg = T 5 + E r (T 2 - T 6 ) 


(C35) 


and 


T 3 



( T 6 - T 5> 


The heat to be transferred is given by 

Q cc mc p (T 2 - T 3 ) 

and by 

Q cc UA(T 3 - T 5 ) 

Ratioing and combining yielded 

Q _ mc p< T 2 - T 3> _ U < T 3 ~ T S> 
Q* m*c*(T* - T*) U*(T| - T*) 

But, from equations (C35) and (C36), 



T 3" T 5 1 ~ E r 


(C36) 


(C37) 


(C38) 


(C39) 


(C40) 
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Therefore, 


E. 


E. 



(C41) 


From reference 13, when the small terms due to wall conduction and fouling factors 
were neglected, 


1 1 1 



Assuming 7j g h = rj g c = 1 and h h s h £ yields 



Therefore, 


(C42) 


(C43) 


U cc h (C44) 

The expression for the local-heat-transfer coefficient may be written from the defi- 
nition of the Colburn factor j as 


h = jGc p (Pr)- 2/3 

where for constant y and fixed geometry, 

j (Re) n (ref. 13) 
G oc m 


Re a — 




(C45) 


(C46) 

(C47) 

(C48) 
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i 


(C49) 


c cc 
P 


M_. 


w 


and 


Pr cc ..... M 

kM 

w 


(C50) 


Therefore, 


JL U /m\“ + VM*\“ + 2 / 3 /k\ 2 / f M w \ 1/3 

h*~u*"w*j u y yj 

Substitution of equation (C51) into equation (C41) yielded 

E r _ /^\^f /3 /k\V 3 / tf *\^/y E ? \ 

1-E r ‘ m *'\ M w/ W L/ \1-Ejj 


(C51) 


(C52) 


Recuperator design predictions were used in equation (C52) to solve for n. A value 
of n = -0. 48 was found to fit the predictions. The predicted changes were verified 
during the recuperator air acceptance tests. However, a value E* = 0. 941 was pro- 
jected rather than the design value of 0. 95 given in table I. When E* = 0. 941 was used, 
the expression became 


/M Y°- 48 /m \ 2 / 3 / 

r = 15. 95/ h ' ' w 




A 


1 - E 


1.28, 


‘ 83,8 ^ \7.22X10" 5 


3.2 5x10' 
M 


,0. 19 


E 


1 - E. 


= 15.95 


M,_ V°- 48 /m...\ 2/3 / k \2/34 0 ^„-5N 0 - 19 


> (C53) 


w 


\0. 581 


183.8/ \0. 0450; 


'484x10 


J 


where the reference or starred conditions (eq. (C52)) were evaluated at the flow condi- 
tions for 10 kilowatts of gross power output. The gas properties were evaluated at 


T = 


T 2 + T 5 
2 


(C54) 
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Waste heat exchanger . - A slightly different approach was taken for the waste heat 
exchanger. Its heat-transfer effectiveness is a function of both the gas and liquid mass- 
flow rates, the specific heats, and the heat to be transferred. Expressed functionally, 


E 


w 


f'(C , N , ) 
v w’ t u' 


(C55) 


where 



(C 56) 


and 


N tu S 


UA 


UA 


(C57) 


mC p) m “S 


After E was determined, the liquid-side temperatures resulted from 
w 


T 4, r T 3- 


T 3~ T 4 


E 


(C58) 


w 


T 3 , l " T 4 ,l + C w (T 3 " T 4 ) 


(C59) 


The local gas- side heat-transfer coefficient was assumed to vary as in equation (C51) 
for the recuperator. Based on the design geometry and conditions, h* = 0. 00556 Btu per 
second per square foot per °R (^114 J/(sec)(m 2 )(K)) was calculated. From the design 
predictions, n = -0. 65 was found to fit. Hence, 


h = 0. 00556 
g \1.28 


m \°- 35 /2.01 x10" 5 


, 0 . 02 ; 


83. 8' 


.1/3. 


v2/3 




M 


ML 


w 


U. 77x10 


-6 


h_ = 114 


m \ 0,35 /3. 00 x10 


0. 581 


■5\°- 02 / 83 8 \ 1/3 / i , \2/3 


M 


M. 


w 


0. 0298 


r (ceo) 
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where the reference or starred conditions (eq. (C52)) were evaluated at the flow condi- 
tions for 10 kilowatts of gross power output. The gas properties were evaluated at 


T 3 + T 4 


(C61) 


The local-heat-transfer coefficient on the liquid side may be expressed by equa- 
tion (C45), but in this case 


Therefore, 


i - [(Re),]" 

(C62) 

m ? 

( Re ) z oc _ £ 

H 

(C63) 

G z am z 

(C64) 

C p,Z = f '^) 

(C65) 

II 

'o - ' 

(C66) 



(C67) 


The exponent n was assumed to be the same for both surfaces, and for the design con- 
ditions, 


0.0137 


Btu 

(sec)(ft 2 )(°R) 


= 281 " 

(sec)(m 2 )(K) 
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was calculated, 
tied to 


Changes in and 


were neglected, and the expression was simpli- 


h z = 0.0137 


h, = 2811 


( m l c y,V 

0.0870; 


to,c 7 > 
l P.Z 

165 , 


,0.35 


A 


(C68) 


where the liquid capacity rate was evaluated at design C, . and c 7 is evaluated at 

W p, 6 


- _ T 3, l + T 4 ,1 

1 i 


(C69) 


After h, and h are calculated, the procedure to find E was that of Kays and 
London (ref. 13). The crossflow equations for "unmixed" liquid flow and "mixed" gas 
flow were used. For completeness, the equations, with substitutions for the design 
geometry, are presented. The overall gas-side heat-transfer coefficient results from 


where 


with 


1 - 1 + 1 


U g , s>g h g 0-336 


(CIO) 


1s,g s l- 0-682(1 -H f , g ) 


(C71) 






tan h^0. 268 ft) 

0. 268 |/i T 

tan h(o. 00337 |/lT) 
0.00337 J h 


A 


S 


> 


(C72) 
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and where 


\i= 1 -°- aw-iu) 


with 






tan h^O. 076 |/h^) 
0. 076 |/h^ 

tan h(o. 000957 A) 


0.000957 




The number of thermal units was then evaluated from 






N tu = 


<pmY c 


mc n 


where 


(UA) „ 

c p = 0.93 

(UA)* 


based on the air acceptance test results; and from the design, = 405 square feet 
(37.6 m ). Since there are eight liquid passes, 


N, 


N. 


tu, p 


tu 


Then, 


r p -_i-e' C -Vp 


-T/C 

E p = 1 - e P w 


(C73) 

(C74) 

(C75) 

(C76) 

(C77) 

(C78) 

(C79) 
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1 - C E 
5- w P 

‘- e p 

(C80) 

and 


j» 8 _ I 

E w = , ( C w * !) 

? - C w 

(C81) 

If C = 1, by L'Hospital's rule, 
w 


E . 8E P 
W 1+7E P 

(C82) 

Gas Pressure Drops 


For each heat exchanger, the change in total-pressure drop was assumed to vary as 
the change in core friction pressure drop. The following proportionalities lead to the 
relation in equation (C89): 

Ap ^ f pv 2 
P P 

(C83) 

where 


f (Re) n cc ^ 

(C84) 

and 


V oc — 

P 

(C85) 
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Therefore, 



But 


P 


M p 
w v 


and 


Ap m n T 


M n M w V P 


Hence, the general expression used for each heat exchanger was 


(C86) 


(C87) 


(C88) 


/ Ap\ _ / m \ n+ ^/ u *\ n A^w\ /T \ / p*\2 /a p\* 
, P / \ m 7 \ M / \M w y\T*/\p / [ p / 


(C89) 


All values of (A p/p)* used in the analysis were based on the total design values. 
Air measurements have shown smaller values in some cases. However, for conserva- 
tism, the design numbers were used. Viscosity terms were evaluated at the average 
temperature through the component. However, temperature and pressure were evalu- 
ated, conservatively, at the inlet conditions of the component. 

Recuperator . - A value of n = -0. 60 was found to fit the design predictions. All 
reference or starred values (eq. (C89)) were evaluated at the flow conditions estimated 
to give 10 kilowatts of gross output, in the following equations: 
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'Ap\ 

p L 


m 


,1. 40/ 


0.011 


1.267, 


3. 20x10* 


-0. 60/ T \ 
Is JL F“ % 




0.011 





> 


(C91) 


Waste heat exchanger - The same value, n = -0. 60, was found to fit the design 
predictions. All reference or starred values (eq. (C89)) were evaluated at the flow con- 
ditions for 10 kilowatts gross output in the following equations: 



Source heat exchanger . - For this heat exchanger, a value of n = -0. 20 was used 
in its design analysis. The same value was used in this analysis. The reference or 
starred values were estimated for the 10-kilowatt-output conditions and were rounded 
off. The expression, as used in the program, was 



(C93) 


Ducting. - The friction pressure drop for ducting may be expressed as 


AP K eX 
P 2gp 


(C94) 
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or, as used in the program 


Ap 

P 


Ap = 

P 


K f { 1 545 \/T \/G \ 

eq \M W JUAp) 


V 


8317\/t\/g) 


k M w /\ 2g/\p, 


> 


J 


The following friction factors were used: 

f = — for Re < 6000 (laminar flow) 
Re 


(C95) 


(C96) 


and 


f = 


0.316 


(Re) 


0.25 


for Re > 6000 


(C97) 


which is the Blasius formula for turbulent flow (ref. 14). 

Values for K eq for each segment of the ducting were calculated from equations (C95) 
for the design pressure drops allowed at the 10-kilowatt-output conditions. These 
values are shown in table II. 


TABLE II. - VALUES OF K 

eq 


Ducting from- 

Equivalent duct 
length-to-diameter ratio, 

K eq’ 

dimensionless 

Turbine to recuperator 

47.2 

Waste heat exchanger to compressor 

38. 8 

Compressor to recuperator 

42. 0 

Recuperator to source heat exchanger 

136.0 

Source heat exchanger to turbine 

87. 5 
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Gas Inventory 


For the calculation of total inventory, the gas loop was divided into nine segments. 
The inventory was calculated as 



144 M w 
1545 





> 




(C98) 


The segments and their estimated gas volumes are shown in table in. 


TABLE m. - ESTIMATED VOLUMES OF 


GAS- LOOP SEGMENTS 


Segment 

Volume, V 

ft 3 

m 3 

Duct from turbine to recuperator 

0.326 

0.00923 

Recuperator, hot side 

1.98 

.0560 

Waste heat exchanger 

.714 

.0202 

Duct from waste heat exchanger to compressor 

. 104 

. 00294 

Duct from compressor to recuperator 

. 123 

. 00348 

Recuperator, cold side 

1.27 

.0359 

Duct from recuperator to source heat exchanger 

.184 

.00521 

Source heat exchanger 

. 380 

.0108 

Duct from source heat exchanger to turbine 

.370 

.0105 

Total 

5. 451 

0.1543 


System Electrical Needs 

Part of the gross alternator power output P n is used to control the power system 

u 

itself. Hence, the gross system output power is given by 


p = p _ p 

*sy, G a, G e 


(C99) 


T 
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Figure 12. - Schematic diagram of estimated system power needs. 


The estimated system power needs are shown in figure 12. Hence, the expression for 
control power is 


P = 1. 1625 + 0. 0250 v n (C100) 

G 3,y \JX 

All of the electrical component power needs are based on approximations. The power 
required by the pump-motor assembly is a function of liquid mass-flow rate. Measured 
acceptance test results show this motor power need to vary between 330 to 430 watts. 

The use of 400 watts would correspond to a total coolant-flow rate of 0. 35 pound mass 
per second (0. 16 kg/sec). 


Radiator Areas 

In the analysis of parametric system operation, approximate values of prime radia- 
tor area are calculated. Over the range of pump design mass-flow rates, the calculated 
pump- motor coolant-temperature rise varied from about 0. 5 to 1.5 R°. For conserva- 
tism, a constant temperature rise of 2 R° (1 K) was assumed. Then, neglecting any 
temperature drop between the liquid and the radiator surface, 

T R,i = T 3 ,1 < C101 > 
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and 


T R, o = T 4, Z ~ 2 


T - T - 1 
l R, o ~ l A,l 1 


(Cl 02) 


The prime, or minimum area, need was then calculated from 



ln < t r, j - t s>< t r, 0 + T s> , 


4 " £T s L < t r, o - T S>( T R, i + T s> 


tan" 


■irklV tan- 1 /^-^ 


(Cl 03) 


This expression is a simplified form of the general expression presented in reference 15. 
The missing term, which adds the effect of the coolant heat-transfer coefficient, has 
been assumed to be small. Mathematically, this is true only for an infinite heat-transfer 
coefficient. 
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